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Model predictive control (MPC) has recently been ranked as one of the best 
modern technics of control for high voltage direct current (HVDC) based on 
voltage source converter (VSC) transmission systems, since it is a simple 
controller, understandable and can be implemented easily, also has a 
dynamic and fast response. The classic double closed loop PI controller has 
many disadvantages such as the difficulty in determining the Kp, Ki, and Kd 
parameters as well as the couplings; It is also less robust against systems 
with slow dynamics and/or having a complex structure. This paper exposes a 
model predictive control for a two-level converter. Direct power control on 
rectifier-side based on model prediction and on the inverter-side AC current 
control strategy is proposed. The controller utilizes the identified model to 
anticipate the comportment of the output voltage of every switching state at 
each sampling frequency. At the next sampling interval, the optimum 
switching state is applied using the cost function as a reference point. The 


MPC method was verified by the tools of the MATLAB/Simulink software. 
The simulation results validate the feasibility and effectiveness of the 
presented control strategy. 
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1. INTRODUCTION 

The direct current (DC) for long-distance power transmission began in 1954 once ASEA, build the 
first link with high-voltage direct-current (HVDC) lines. However, the DC’s history appears before that. 
Thomas Edison failed to keep on against his rival the alternating current AC in the 1890s it became clear that 
ac was more efficient at transmitting electricity over long distances. On the other hand, due to the increase of 
electricity consumption through the world, the DC is nesting again between the branches granting 
a new spot [1]. Nowadays, power electronics are used to convert and control electrical energy. The growing 
consumption of electrical energy has prompted the development of new standards for evaluating energy 
quality and efficiency [2], [3]. Faced with the novelty of semiconductors, a whole new generation of control 
diagrams has been traced [4]. In the 1990s, the growth of semiconductors such as insulated-gate bipolar 
transistor and gate turn-off thyristor, had reached values allowing them to be used by voltage source 
converters (VSC). The very early VSC-based HVDC transmission line with a 50 MW underground cable was 
operational in 1999 on the island of Gotland; there are several names for these transmission cables such as 
HVDC plus, VSC-HVDC, and HVDC light [5]-[7]. The HVDC system requires the elaboration of rigorous 
control techniques of both active and reactive power flows, because it has several advantages such as a 
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decoupled and fast control of the power flow, the improvement of the power quality indicator and supply of 
isolated loads [8]-[13]. 

The model predictive control (MPC) shows a huge advancement in modern automatic control 
theory. At first, the process industry was the perfect field to study and apply the MPC, where it has been in 
use for many years ago [6]. Today, other fields put a lot of attention against the predictive control, such as 
electrical drives and power electronics [14]-[18]. The growing use of MPC is due to the better mathematical 
models for predicting the comportment of the controlled variables in electromechanical systems [14]. 
Moreover, the recent cutting edge technologies microprocessors can perform at high speed and low cost a 
large number of calculations required in MPC [19], [20]. The MPC method has been used in several power 
converters for its high dynamic performance, simple design and ability to regulate the AC signal directly 
without frequency conversion [21]-[26]. 

Recently, some recent research has found very interesting results of MPC predictive control. The 
article [18] draws up a table of some differences between the model predictive controller and the 
proportional-integral regulator and demonstrates the robustness of the predictive control through theoretical 
and experimental calculations. Zhang et al. [27] gives the mathematical modeling and the design of the MPC 
control strategy. Zhao et al. [28] details the MPC control of VSC-HVDC that powers passive networks. 

In this article, we started with the discrete mathematical modeling of the two sides of the VC-HVDC 
system. Direct power control was applied on the rectifier part, while the alternative control based on a 
prediction model was used on the inverter part. In its simplest structure, the MPC control is more efficient 
due to the anticipation of the next values of the variables, which makes it more robust and more stable than 
the classical double closed-loop proportional integral (PI) control. The dynamic response speed is rapid, and 
the steady state and dynamic performance are both good. Finally, based on MATLAB/Simulink platform, we 
have built a VSC-HVDC system that powers a passive network, in order to confirm the suggested control 
strategy. 


2. MATHEMATICAL MODEL OF THE VSC-based HVDC SYSTEM 

Figure | gives the general structure of the VSC-based HVDC system. It consists mainly with three 
main parts: sending end, transmission line and receiving end. Where Table 1 shows in detail the components 
that consists in that system. 


Table 1. System symbols designations 


Item symbol Meaning 

Zyj Equivalent resistance of the AC system source 

L; Commutation reactor inductance 

Ri Commutation reactor resistance 

Uj Three phase AC voltage 
ij Three phase AC current 

ujde DC voltage 

lde DC current that passes the DC link transmission line 
C DC capacitor 


Where j=1, 2 (1: rectifier station, 2: inverter station) 


Sending end Receiving end 


Figure 1. The general structure of VSC-HVDC system 


2.1. Mathematical model of the rectifier side 
According to Figure 1, (1) is determined by the KVL law, in a stationary three-phase coordinate 
system, a continuous mathematical model of the rectifier side can be obtained. 


diik _ š 
Lı ie “kT Rilik — Urk (1) 
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Where k=a, b, c. To decouple currents and voltages, we change the frame by performing Clark transformation. 


di, 


1e “17 Rii — u, (2) 


i u 
Where i = | , u = | a 
1 


u 
ú | and u, = el The current i; can be expressed as a function of the 
1B rp 


current i, as: 

i = 2/3 lia +a ip + a? ie) = ha t firg (3) 
where a = e/"/3). Also, the voltages u, and u,are defined as (4). 

uy = 2/3 (Uia + a Urp + a? uic) (4) 
The voltage u,is given by Utdc as: 

Ur = Srec * Urde (5) 


where Spec is the switching state vector of the rectifier expressed by (6). 
Srec = 2/3 (Sia Hasip t a? Sic) (6) 


In this formula, Six is the switching function of the rectifier in the k-phase bridge arm, which is designed to 
work as mentioned in Table 2. 
Discretization of (2): 


RıTs 
Li 


ilk +1) = (1-45) + Bf) - u] (7) 


The discretization is obtained by approximating the derivative as the difference over one sampling period (T,): 


di _ i(k+1)-i(k) 


x ID- (8) 


dt Ts 


Table 2. Rectifier switching function 
Sık_  Upperarm Lower arm 
0 OFF ON 
1 ON OFF 


2.2. Inverter side mathematical model 
When the AC system is stable (Figure 1), the mathematical model of the inverter side can be 
obtained as: 


dizk _ r 
2 Mak T Relax + Uik (9) 


where, the output voltage of the inverter is given by: 


1 
Uik = U2dc (Soe es nsabe Som) (10) 
By performing the inverse Clark transformation, the af coordinate system of (9) can be obtained as: 


2 = = —Uz — Ratz + UgacSiny (11) 
where Siny = 2/3 (S2a + @ Sz, + a? S2¢). Discretization of (11) using the first-order forward difference 
method yields: 

R2Ts 
L2 


i(k +1) = (1 - #5) ik) + P fuk) - u] (12) 
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3. CONTROLLER DESIGN 
3.1. Rectifier side controller 
Instantaneous predicted input powers are calculated by: 


Pin (k + 1) = Refu, (k + G(R +1)} = talia + gig (13) 
Qin(k + 1) = Imfuk + DG + 1)} = giia — talig (14) 


where i, (k + 1) is the anticipated input current. About the grid fundamental frequency of the network, for a 
small sampling frequency, we can assume that u,(k +1) =~ u,(k). A proportional integral regulator 
controls the intermediate circuit voltage. The control signal of the controller symbolizes the power required 
to balance the voltage error of the intermediate circuit. 


* ki * 
Pin = (kp + Uae > Uzac) + Piac (15) 


Where Piac = UWraclac + Pic, here Picis the active power demanded by the capacitor. At steady state we 
suppose that Pje ~ O thus Piac = uUygcige. So as to stabilize the powers at their setpoint values, the 
following objective function is expressed by: 


Irec = rn ane 1)| + |Ph — Pin(k + 1)| (16) 


where Ph and Q;, are the active and reactive power setpoints respectively. The cost function grec 
summarizes the rectifier’s desired performance: control the active power Pin and reduce the reactive power 
error. The block diagram and the designed Simulink model of the control strategy are given in Figures 2(a) 
and 2(b) respectively. 


MPC Controller for VSCI (Rectifier) 


(b) 


Figure 2. Predictive power control scheme for the rectifier: (a) block diagram and (b) Simulink model 
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3.2. Design scheme of rectifier side controller 


The purpose of the current control is to reduce the static error. To satisfy this condition, (17) can be 
wrote as a cost function, which calculates the error between the references and the anticipated currents: 


Jinv = lita(k + 1) ~ ifa (k + 1)| + 


iig(k + 1) — Pplk +1) (17) 


where if,,(k +1) and i7 p(k + 1) are the real and imaginary parts of the predicted load current vector 


iP (k + 1). By assuming that i*(k + 1) ~ i* (k), this supposition can give a one-sample delay in the reference 
pursuit, which is not a problem if a high sampling period is considered. The block diagram and the designed 
Simulink model of the current control strategy for the inverter is given in Figures 3(a) and 3(b) respectively. 

The basic steps of the predictive current control algorithm are as follows, see Figure 4: i) measure 
the controlled current (i2(k)) at the kth sampling period; ii) carry out the best switching state considered in the 
previous sampling frequency; iii) predict the output currents for all possible switching states for the 
upcoming (k+1)th sampling instant, based on the discrete time model; iv) determine the cost function for the 
predicted output currents based on the control objectives; and v) choose the switching state that decreases the 
cost function and stock it for using at the next sampling period. 


i(k) 


i(k+1) ae 
Of inv 


g function 
i(k) Predictive => 


model 


<signal4> 


<signal6> 


MPC Controller for VSCII (Inverter) 


(b) 


Figure 3. Predictive current control: (a) block diagram and (b) Simulink model 
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Measure i (k), in(k) 


Predict 
iP (k+) if (k+1)2 
s---- iP (k+1)7 


Evaluate 
gir (k+1):); gG? (k+) 


Figure 4. The predictive current control algorithm 


4. SIMULATION AND RESULTS 

The simulation model of VSC based HVDC system is built in MATLAB/Simulink (with the help of 
MATLAB embedded system functions) to validate the model with predictive control method. The schematic 
layout used for simulation is given in Figure 5. The first converter (VSC-I) used to regulate powers, whereas 
the second converter (VSC-II) used to trace the AC current. Between two converters, a 100 km long 
transmission line model is also used; the simulation model is also built under the following parameters: The 
three-phase sending end line (base) voltage is set to 230 kV, rated (base) voltage of DC-side widce=360 kV, 
DC capacitor Cae=6000 uF. The receiving end line (base) voltage is set to 230 kV; the base apparent power 
800 MVA, the system frequency is set to 50 Hz and the sampling period is set to 25 usec. 


Terminal-I 


[occececcecceceeee Terminal-Il 

' : pee ee 

H ' 

i Co | 

Por E “ H 
AC-Grid Bmw ewe wenn enseee 


VSC-I Line VSC-II 


Figure 5. VSC-HVDC system used in our simulation 


4.1. Normal operation 

The AC voltage of rectifier side is fixed to 1.0 pu (230 kV). The DC voltage is fixed to 1.0 pu 
(360 kV). The reactive power reference value is set to 0 pu (0 Mvar) and the active power reference value is 
set to 0.75 pu (600 MW). The following simulation waveforms are acquired at steady state. Allowing to 
Figures 6 and 7, the system in steady state has the following features: 
—  Atrectifier side 

The control approach of the rectifier part is to keep DC link voltage, active and reactive power on 
their setpoints values. It is identified from Figure 6, the first 0.1 s is the commencement process, and the DC 
link voltage slowly becomes stable at 1.0 pu that be given as shown in Figure 6(c). Also, Figure 6(d) illustrates 
that the active and reactive power are also stable at their given values. In addition, Figures 6(a) and 6(b) shows 
the input voltage and current respectively, in which they are sinusoidal waveforms and are also in phase. 
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Pin and Qin (pu) 
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Time (sec) 


(d) 


Figure 6. Rectifier side: (a) AC input voltage, (b) AC input current, (c) DC link voltage, and (d) input active 
and reactive power 
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— At the inverter side 

The AC output voltage is sinusoidal signal and in phase with the output current iz. The output 
current iz tracks its reference value (wich is around 0.74 pu, since we take in consideration the transmission 
line losses), as shown in Figure 7(a). Also, Figure 7(b) demonstrates that the active power is slightly less than 
the rated value; this is obvious because line parameters can influence the transmitted power over the DC link 
(since we are transmitting power to a passive network i.e there is no compensation for the power losses). The 
output voltage u2 is almost a sinusoidal waveform and in phase to the output current as mentioned in 
Figure 7(c). 


4.2. Malfunction operation 

In order to test the dynamicity and the fast response of the designed controllers, a disturbances/step 
response must be applied to a part of the whole system. The goal from those tests is to see the behavior of the 
fault’s applied section (eg: the sending or receiving end) and its impact to the other end and vice versa. 
According to the previous arguments, we apply two step tests, the first is a step in the input AC voltage of 
-0.2 pu at t=1 s; and the second one is a step of +0.25 pu in the reference output current at t=1 s. 


Miesa 


Teens 
(a) 


I2 andi2ref (pu) 


Time (sec) 


(b) 


u2 (pu) 


0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37 0.38 0.39 04 
Time (sec) 


(c) 


Figure 7. Inverter side: (a) AC output current, (b) output active power, and (c) AC output voltage 
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4.2.1. Step in the input AC voltage (u1) 

After we apply a step of -0.2 pu in the input voltage (u1) at t=1 s (the fault is cleared after three 
cycles), the following results are obtained. Referring to Figure 8 we can deduce the following results: 
—  Atrectifier side. 
After the fault is applied at t=1 s, the input voltage and current present a fast dynamic behavior to keep the 
controlled variables (the active and reactive power) at their reference trajectories seen from Figures 8(a) and 
8(b) respectively, the active power showed two peaks of around -0.15 and +0.15 pu at the start and the end of 
the applied fault respectively, after the fault is cleared the active power restored to its stable state at 0.75 pu 
as shown in Figure 8(c), while the currents are kept sinusoidal and in phase with the supply voltages, even 
during the transient. 
- At the inverter side. 
All the waveforms at the inverter side is kept unchanged like the ones shown in Figure 7 (since the power 
transmitted over the DC link is kept almost constant during the fault). 
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Figure 8. Comportment of the predictive power control in the input AC voltage of —0.2 pu at t=1 s (a) AC 
input voltage, (b) AC input current, and (c) input active and reactive power 
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4.2.2. Step in the reference output current (i*2) 

After we apply a step of +0.25 pu in the reference output current at t=1 s, the following results are 
obtained. The comportment of the predictive control scheme of the load current reference is given in 
Figure 9. The inverter currents and active power react with fast dynamics to this variation (shown in Figures 
9(a) and 9(b) respectively) that leads the rectifier to compensate the power needed, this is done by modifying 
the input active power setpoint. In addition, Figure 9(c) shows that the input currents has a fast dynamic 
performance. While Figure 9(d) demonstrates that the input active power has a good performance due to this 
variation while the transient does not disturb the reactive power. 


i2 and i2ref (pu) 
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Figure 9. Comportment of the predictive power control (a) output AC current, (b) output active power, 
(c) AC input current, and (d) input active and reactive power 


5. CONCLUSION 
In this article, we have discussed the discrete mathematical model of the inverter and converter parts 
of the VSC-HVDC system; as well as their predictive control that is used to control and maintain the 
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controlled variable in their trajectories. Knowing that the predictive control uses a finite set of switching 
states: i) that can be produced by the power converters; ii) that are employed to anticipate the future 
comportment of the controlled variables (in our case the direct power flow at the rectifying station and a AC 
current at the inversion station) for each switching state; iii) the simplicity of the control structure since there 
is no complicate calculation or algorithms, also there is no inner loop current control which improves the 
system reliability; iv) fast dynamic response especially at faults, which gives the proposed control method a 
good advantage rather than the classical methods used before (eg: PI control); and v) good steady state 
performance. 
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